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I. 


Introduction 


The  use  of  a  coherent  radar  for  tracking  a  moving  object  has  long  been 
established  as  an  important  technique.  By  measuring  the  Doppler-phase,  the 
velocity  and  acceleration  of  a  target  may  be  determined  according  to  some  pro¬ 
cedure  such  as  matched-filter  estimation.  Such  estimations  usually  contain  the 
implicit  assumption  that  the  target  is  a  point  scatterer.  However,  most  targets 
are  extended  and  exhibit  complex  patterns  of  scattering  amplitude  and  phase  as 
a  function  of  viewing  (aspect)  angle.  Thus,  questions  often  arise  as  to  the  effect 
of  the  scattering  phase  of  a  tumbling  or  oscillating  target  on  the  estimations  of 
velocity  and  acceleration. 

In  discussing  the  effects  of  scattering  phase  on  velocity  and  acceleration 
estimates  there  are  two  possible  viewpoints.  One  might  be  performing  an  ex¬ 
periment  where  specific  knowledge  of  the  scattering  phase  versus  target  orien¬ 
tation  was  pre-measured  and  target  orientation  versus  time  is  determined.  In 
this  case,  one  can  subtract  the  scattering-phase  effect  and  thus  obtain  better 
trajectory  estimates.  On  the  other  hand,  one  may  simply  wish  to  establish  a 
reasonable  bound  for  the  effect.  It  is  then  only  necessary  to  combine  an  estimate 
of  the  rate  of  tumbling  or  oscillation  with  some  knowledge  of  the  scattering 
phase.  It  is  the  second  viewpoint  which  we  shall  pursue. 

In  Section  II  we  establish  a  specialized  model  for  the  phase  of  the  re¬ 
ceived  signal.  In  Section  III  we  briefly  discuss  the  estimation  procedure  for 
a  point  target.  By  assuming  the  angular  rotation  of  the  target  is  sufficiently 
small  during  the  measurement,  we  may  show  that  the  scattering  phase  intro¬ 
duces  bias  in  the  estimates  of  velocity  and  acceleration.  The  results  are  de¬ 
rived  for  the  procedure  of  fitting  a  parabola  to  the  Doppler  phase.  Finally, 
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in  Section  IV  we  show  an  example  of  measured  scattering-phase  data  and  in¬ 
dicate  the  velocity  error  one  might  encounter  if  the  target  were  not  stabilized. 

II.  Scattering  Phase  Model 

The  phase  of  a  received  signal  is  (for  transmitted  and  received  polar¬ 
izations  i  and  j  respectively) 

^•(t)  =  - 0UQ  [  t  -  A(t)  ]  -  ^j[t  -  A(t)  1  f 

where  A(t)  is  the  round-trip  delay  to  the  target  and 
phase  evaluated  at  the  time  the  signal  was  scattered, 
lated  to  the  position  of  the  target  by 

cA(t)  =  2  |  T[t  -  lift)]  I  .  (2.  2) 

The  time  dependence  of  is  implicit  through  its  dependence  on  target  orien¬ 

tation.  This  orientation  is  described  by  three  Euler  angles  in  the  radar  coordi- 

A 

nate  system  (line-of-sight,  k  ,  and  two  orthogonal  linear  polarization  vectors, 

A  A 

Gj  and  e^).  For  the  special  case  of  an  axially  symmetric  target  only  two 
angles  are  necessary  and  it  is  convenient  to  choose  them  as 

A  A 

cos  0  =  k  •  s  (2.  3a) 

A  A 

sin  0  cos  c/?  =  Gj  *  s  (2.  3b) 


(2.  1) 

is  the  scattering 
The  delay  A(t)  is  re- 
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where  s  is  a  unit  vector  along  the  symmetry  axis.  It  can  further  be  shown 
that  the  scattering  matrix  for  arbitrary  cp  can  be  calculated  from  the  special 
cases  of  cp  =  0  and  cp  -  tt/2  (i.  e.  in  or  orthogonal  to  the  plane  of  in- 

A  A 

cidence  defined  by  k  and  s  ).  We  may  then  specialize  Eq.  (2.  1)  to  the  case 

for  which  5. .  depends  on  time  through  only  one  variable,  0  [t  -  i  A(t)  ]. 
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III.  Estimation  of  Velocity  and  Acceleration 

For  a  point  target  the  Doppler-phase  is  determined  by  the  range  delay 
A(t)  .  Equation  (2.  2)  indicates  that  A(t)  contains  a  measure  of  the  target 
trajectory;  in  particular,  the  radial  velocity  and  acceleration  are  to  be  esti¬ 
mated.  In  order  to  simplify  the  analysis  let  us  consider  the  special  case  for 
which  the  trajectory  is  coincident  with  the  radar  line  of  sight.  (In  the  general 
case  it  is  necessary  to  deal  with  radar  coordinates  including  angles  and  angle 
rates .  ) 

If  the  scattered  signal  is  first  received  at  t  =  T  ,  we  expand  A(t) 
and  r(t)  in  series  about  this  time.  For  the  special  trajectory  being  considered, 
these  expansions  are  probably  valid  for  a  long  interval  of  time;  however,  for 
non-zero  angle  rates,  there  may  be  a  more  limited  interval  of  convergence 
We  assume  that  the  series  are  certainly  valid  over  the  interval  of  time  nec¬ 
essary  for  any  specific  method  of  data  processing  (typically  on  the  order  of 
.  05  to  .  5  second  for  good  velocity  resolution).  Assuming  the  expansion 

A(t)  =  T  -  B(t  -  T)  -a(t  -  T)1 2  +  .  .  .  ,  (3.  1) 

1.  S.  L.  Borison,  "Diagonal  Representation  of  the  Radar  Scattering  Matrix 
for  an  Axially  Symmetric  Body",  IEEE  Trans,  on  Ant.  and  Prop.  ,  Vol  AP-13, 

No.  1,  p.  17  6,  January  1965  (DDC  613887). 
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we  use  Eq.  (2.  2)  to  relate  T,  3,  and  a  to  range,  radial  velocity  and  acceler- 

2 

ation.  The  results  are 


T 


0 


a 


where 


v 


a 


|  r(T/2) 

2y/c 
1  +  v/c 

_ ih _ 

(1  +  v/c)3 


dr(t)  | 
dt 

d2r(t) 


t  =  T/2 


T/2 


(3.  2a) 


(3.  2b) 


(3.  2c) 


(3.  3a) 


(3.  3b) 


Thus,  estimates  of  T,  3  »  and  OL  are  equivalent  to  estimates  of  r,  v, 
and  a  at  the  time  T/2. 

The  specific  details  of  an  estimation  procedure  will  depend  on  the  dur- 

2 

ation -bandwidth  product,  TB,  of  the  signal.  For  small  TB  it  may  be 
sufficient  to  estimate  the  velocity  and  derive  acceleration  as  the  time  deriva¬ 
tive  of  v  ;  for  large  TB  ,  it  may  be  necessary  to  simultaneously  estimate 

A  A 

v  and  a  .  Whatever  the  scheme  we  shall  find  the  estimates  8  and  a 
with  tolerances  (or  errors)  6  8  and  6a  respectively.  (Neglecting  terms  in 

2.  E.  J.  Kelly  and  R.  P.  Wishner,  "Matched-Filter  Theory  for  High-Velocity, 
Accelerating  Targets",  IEEE  Trans,  on  Military  Electronics,  Vol  MIL-9, 

No.  1,  pp.  56  -  69,  January  19  65.  (DDC  616440). 
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A  A  * 

/c  c 

c,  Eq.  (3.  2b)  implies  that  v=  —  (3  and  6v  =  -r-60).  We  now  wish  to  determine 

2  ^ 

the  additional  errors  incurred  due  to  phase-center  motion. 

The  method  for  calculating  the  effects  of  phase -center  variations  will  de¬ 
pend  on  the  duration-frequency  product,  TF  [  F  is  essentially  the  instantaneous 
value  of  d0(t)/dt].  For  TF  »  1,  the  phase,  E  [  0  (t)  ]  ,  will  undergo  many  re- 

A  /\ 

petitive  cycles;  calculation  of  the  errors  induced  in  v  and  a  will  tend  toward 
zero,  independent  of  the  estimation  scheme  due  to  long-time  averaging.  For 

A  A 

TF  ~  1,  the  errors  induced  in  v  and  a  will  depend  on  the  initial  value  of  £ 
and  on  the  particular  method  of  data  processing.  For  TF  1  the  errors 

A  A 

induced  in  v  and  a  depend  only  on  the  local  values  of  ^  and  its  time  deriv¬ 
atives;  furthermore,  these  errors  may  be  determined  for  the  procedure  of 
fitting  a  parabola  to  the  Doppler  phase. 

The  method  to  be  followed  is  straightforward.  For  TF  «1,  we  assume 
5  may  also  be  expanded  in  a  series  about  t  =  T.  The  terms  in  (t  -  T)n  may  be 
associated  with  corresponding  terms  from  the  expansion  of  A  (t) .  It  is  obvious 


that  the  estimates  of  8  and  a  are  simply  biased  by  the  changes  in  E;  these 
biases  are  independent  of  the  estimation  procedure,  and  the  error  estimates 

A  A 

63  and  6a  are  unchanged.  The  results  of  expanding  5  in  a  series  (neglecting 
terms  of  order  v/c  )  is 

§  [  e  {  t  -  4  A  (t)  }  ]  =  §  [  9(T/2)  ]+  (t  -  T)  —  —  l^/o  + 


de 

de  9<T/2>  dt 


T/2 


1  ,  .2  rd^F  |  /de  %  2  d? 

7  *'T  [^-2  1  a  (T/2)  <—  > 


(3.4) 


dt 


(T/2) 


(d  o  L  a  de.  ]+. . 


de 


S(T/2)  '^1  +  -  — ) 


dt" 


c  dt 


T/2 


*  Hereafter  we  shall  drop  the  polarization  subscripts. 
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By  associating  terms  from  Eqs.  (3.  1)  and  (3.4)  in  Eq.  (2.  1),  it  is  clear  that 

A 

the  8  estimate  will  be 


:  _  2v  1  d?  |  d0  | 

=  C  „  de  6{t/2>  It/2 

o 

Thus  the  Has  in  the  velocity  estimate  is 


bias 


_ c _  d£ 

2uu  #  d0 


0(T/2) 


de 

dt 


(3.5) 


Similarly  Eqs.  (3.  1)  and  (3.4)  lead  to  an  evaluation  of  the  instantaneous 
acceleration  bias. 

To  use  Eq.  (3.  5)  it  is  necessary  to  have  some  knowledge  of  the  scatter- 
ing  phase  and  of  the  body  dynamics.  If  we  assume  a  reasonable  value  for 
d@/dt  =  2TrF(t)  ,  where  F(t)  has  only  small  variations,  we  may  calculate 


Vbias 


\F  d£ 
2  d0 


for  any  aspect  angle  ,  0.  For  the  acceleration  bias  we  find 


(3.  6) 


=  ir\FZ  [d2?/de2  +  — -  (|+^T)d?/de]  .  (3.7) 

2ttF  dt 

(Note  that  akias  depends  on  the  true  value  of  a.)  In  some  cases  ?(0) 
may  be  available  by  static -range  measurements;  if  not,  some  reasonable 
theoretical  predictions  can  be  made. 
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As  a  final  point  in  this  section,  let  us  note  that  a  consistency  check  is 
available  for  using  Eqs.  (3.  6)  and  (3.  7).  For  the  observation  (integration) 
time  t  ,  A0  ^  tF  .  If  the  function  ?(0)  can  be  well  represented  by  only 
linear  and  quadratic  terms  over  any  interval  of  length  A  0 »  then 
Eqs.  (3.  6)  and  (3.  7)  are  valid.  If  A0  (i.  e.  TF)  just  encompasses  a  region 


of  rapid  linear  (or  quadratic)  change  in  £  ,  the  estimate  of  v^ag  (or 

A 

^ias  ^  nearly  maximized.  If  A0  is  too  large  (i.  e.  tF  too  large) 

Eqs.  (3.  6)  and  (3.  7)  will  tend  to  overestimate  by  assuming  large  local  rates 
of  change  in  5(8)  • 


IV.  An  Example 

As  an  example  of  data  for  use  in  Eqs.  (3.  6)  and  (3.  7),  Fig  1  shows  the 

scattering  phase  measured  at  L-band  for  a  metal  cylinder  (37.  345  inches 

3 

long  and  4.  194  inches  diameter).  Note  that  there  is  no  ambiguity  in  the  phase 
since  one  can  readily  identify  the  fly-back  of  the  pen-recorder.  A  3  60°  change 
in  phase  corresponds  to  a  deflection  from  the  center  to  the  outer  circle  in  Fig 
1.  Figure  2  shows  the  corresponding  radar  cross  section.  Phase  calibration 
data  are  also  presented  in  Ref  3,  however,  the  data  is  sufficiently  linear  for 
our  purposes. 

A  brief  inspection  of  Fig  1  shows  that  quadratic  phase  changes  might 

be  seen  for  intervals  of  data  such  that  0  <  A0  =  TF  10°  .  Except  near 

9  =  0,  it/ 2  ,  tt  and  3tt/2  ,  if  20°  &  A0  <  90°  ,  ?(0)  could  be  fitted  by  a 

dominant  linear  dependence  - though  not  derivable  from  local  behavior 

3.  N.  Landry,  ML-Band  Phase  Measurement  System11  TCR  No.  49,  Radio  Corp. 
of  America,  Missile  and  Surface  Radar  Division,  (16  January  1964).  (U) 
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(i.  e.  dF  /d0).  What  we  learn  is  that  Eqs.  (3.  6)  and  (3.  7)  should  be  valid  for 
A0  ^  10°  for  any  0  .  For  10°  -S  A0  5  20°  ,  estimates  of  velocity  and  accel¬ 

eration  errors  depend  on  the  initial  value  of  0  and  the  details  of  the  data  pro¬ 
cessing.  For  20°^  A0  ^90°  the  specific  value  of  0  is  important 

A 

however,  application  of  Eq.  (3.  6)  will  probably  yield  an  upper  bound  for  v*bias 

because  dF  /d0  will  usually  overestimate  the  linear  behavior  of  F(9)  for 

o  A 

this  interval.  For  A0  »  90  ,  the  phase  will  be  cyclic  and  errors  in  v 

A 

and  a  will  probably  tend  toward  zero  due  to  smoothing. 

To  be  specific  let  us  determine  the  maximum  velocity  bias  as  a  function 
of  TF  for  this  target.  For  tF  ^10°  we  simply  locate  the  maximum  value 
of  dp/d0  .  One  can  see  that  large  values  of  dF  /dB  occur  four  times  for 
each  90°  change  in  0  .  The  maximum  value  of  (A?/A0)  is  on  the  order 
of  50°/l°  =  50.  Thus,  by  using  Eq.  (3.  6)  regardless  of  the  value  of  TF, 
we  would  estimate  (X  ~  1  ft  at  L-band) 

max  |  v^.  I  ^  25  F  ft/sec 

for  F  in  sec  ^  .  In  fact,  for  a  realistic  case  of  T  =  .  05  sec  and 
F  =  .5  sec  ^  ,  TF  <  10°  ,  and  we  might  well  observe  a  velocity  bias  on 
the  order  of  25  ft/sec.  This  value  would  be  considerably  larger  than  the  vel¬ 
ocity  resolution  of  the  signal  (X/2T  10  ft  /sec). 
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Figure  2.  Cylinder  (cross-section)  1300  MCS,  VV  Polarization. 
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